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Abstract 

b 

1 We study the scale dependence of twist-3 distributions defined with chirality-odd quark-gluon 

operators. To derive the scale dependence we explicitly calculate these distributions of multi-parton 
states instead of a hadron. Taking one-loop corrections into account we obtain the leading evolution 
kernel in the most general case. In some special cases the evolutions are simplified. We observe that the 
obtained kernel in general does not get simplified in the large- iV c limit in contrast to the case of those 
twist-3 distributions defined only with chirality-odd quark operators. In the later, the simplification 
O ; is significant. 

> 

\q | Predictions can be made from QCD with the concept of factorizations for processes with large momen- 

tum transfers. A typical example is DIS. For unpolarized DIS, the differential cross section at the leading 
power of the momentum transfer Q is predicted as a convolution of perturbative coefficient functions 
with parton distribution functions(PDF's). PDF's are defined with twist-2 QCD operators and describe 
nonperturbative effects of hadrons. Although PDF's can not be predicted with perturbative QCD, but 
they scale dependence, hence the Q 2 -dependence of the differential cross section, can be determined by 
perturbation theory. The dependence is governed by the famous DGLAP equation. In the past, the 
predicted scale-dependence or DGLAP equation has played an indispensable role for testing QCD as the 
^ . correct theory of strong interaction. 

In general, factorized differential cross sections also contain contributions involving hadronic matrix 
elements of higher-twist operators. Although these contributions are suppressed by inverse powers of Q, 
they contain more informations about inner structure of hadrons. Among them, the most interesting 
are those involving twist-3 operators. These contributions are responsible for certain asymmetries in 
differential cross sections. These asymmetries can be measured in experiment and hence provide infor- 
mation about hadronic matrix elements of twist-3 operators. A well-known example is the study of Single 
transverse- Spin Asymmetries(SSA). The asymmetries can be factorized with the ETQS matrix elements 
defined with chirality-even quark-gluon operators at twist- 3 [U [2]. The scale dependence of these twist- 
3 operators have been studied in O HI El [8] . Besides them, there are chirality-odd quark-gluon 
operators at twist-3. In this work, we study the scale dependence of these operators. 

We consider a spin-1/2 hadron moving in the z-direction with the momentum = (P + , 0, 0). 
We will use the light-cone coordinate system, in which a vector a^ 1 is expressed as = (a + ,a~ ,a±) = 
((a + a 3 )/\/2, (a — a 3 )/\/2, a 1 , a 2 ) and a\ = (a 1 ) 2 + (a 2 ) 2 - In this coordinate system we introduce two 
light-cone vector: = (0, 1,0,0) and = (1,0,0,0). There are two distributions can be defined with 
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chirality-odd quark-gluon operators at twist-3. They are: 



f F (x 1 ,x 2 ) = g s J d ^^ e - i y^ p++i y^ p \P\^{y 1 n)(i 1 ^ 1 + ) G^(0)^(y 2 n)\P) , 

Xf A ( Xl ,x 2 ) = 5s |^^ e -^ lP+ +™ P V^^ (1) 

where the matrix element in the first line is spin-averaged and that in the second line is of a longitudinally 
polarized hadron with the helicity A = ±1/2. x\ y2 are momentum fractions. The definitions are given in 
the light-cone gauge n ■ G = 0. In other gauges gauge links along the direction n should be supplemented 
to make the definitions gauge invariant. From general principle one can show that the two distributions 
are real and: 

f F (x 1 ,x 2 ) = f F (x 2 ,x 1 ), f A (xi,x 2 ) = -T A (x 2 ,xi). (2) 

Replacing the gluon field strength tensor G + ^(0) in Eq.([TJ) with the covariant derivative = + 
ig s G^(0), one can obtain another two twist-3 distributions. With equation of motion one can relate those 
two distributions to the two distributions defined in Eq.Q, respectively (see, e.g., [§1 [TU1 [T4"])- 

The operators in Eq.(pQ) contains the operator of the gluon field strength tensor. There are two 
chirality-odd twist-3 operators which only consist of quark field operators. Correspondingly, one can 
also define two distributions. One is e(x) for unpolarized hadrons, another is h F (x) for longitudinally 
polarized hadrons. Again, these two distributions are not independent. One can use operators identities 
in or equation of motion to show that these two are related with the above two and plus some 
contributions with local operators. The evolution of e(x) and h F {x) have been studied in [I2j [I3j [14] 
and the evolution equations have been solved in moment space. The evolution of twist-3 quark-gluon 
operators has been studied in [15J with the emphasis on the solutions of evolution equations. In this work 
we derive the evolution kernel for the two twist-3 distributions defined in Eq.([T]) in momentum fraction 
space with a different method. 

Under renormalization there is no mixing between the two operators in Eq.([T]). The evolution kernels 
can be conveniently expressed with one function by introducing the combinations: 

f±(xi,x 2 ,n) = {P,±\0(xi,x 2 )\P,±) = f F (xi,x 2 ,fJ,) ±f A (xi,x 2 ,fJ,), 

0( Xl ,x 2 ) = 5s |^^ e -«^^ + +^^ + ^(y 1 n)i 7±M 7 + (l+ 75 )G^(0)V'(y 2 n). (3) 

The functions T±{x\ } x 2 ) are nonzero in the region |xi,2| < 1 and \x\ — x 2 \ < 1. The scale-dependence 
can be written in the form: 

^-^T±(a;i,X2,M) = J d£id£ 2 J r ±(x 1 ,x 2 ,£, 1 ,&)f±(£i,£ 2 ,fi). (4) 

The integration region of £1,2 is fixed by the support of T±. It is easy to find that F + is related to 
Here we will take T- to give our results. 

The distributions are defined for a given hadron, but the kernel does depend on the hadron. It is 
completely determined by the dynamics of QCD. For large [/, it can be calculated with perturbative 
QCD. Because of this one can use various parton states to calculate the distribution T_ to find its \i- 
dependence, hence the kernel For the case with operators of twist- 2 one can use single-parton state 
for the purpose. But for the two distributions defined here, one simply finds null results with single-parton 
states because the chirality is flipped by the operators. Therefore, one has to use multi-parton states 
to calculate the distribution. By using multi-parton states factorizations of SSA with twist-3 operators 



2 



have been studied in [TU [T71 [18] . In [7] such multi-parton states have been employed to study the scale 
dependence of twist-3 operators relevant for SSA. 

We introduce the following state \n(\)) as a superposition of single- and multi-parton states: 

|n(A)) = \q) + Cl \qG) + c 2 \qGG) + c 3 \qqq) + ■■■. (5) 

The state |n) is with the momentum (P + , 0,0,0). All partons in the partonic states moving in the in- 
direction and the sum of their momenta is that of |n). The sum of helicity of partons is A. In the above 
we have suppressed quantum numbers of partons, which will be specified later. In principle one can 
introduce wave functions depending on momenta of partons. For simplicity, we take these wave-functions 
as 5-functions and hence Ci(i = 1,2,3- ••) are real constants. If we calculate T_ of the state |^(A)) 
instead of a hadron, we obtain nonzero contributions from interference between different partonic states. 
At tree-level, the contributions can be schematically written as: 

f_(xi,x 2 ) = C 1 (q(-)\0(x 1 ,x 2 )\q(+)g(-))+C 2 (q(-)g(+)\0(x 1 ,x 2 )\q(+)) 

+C3^(-)g(-)|0(x 1 ,x 2 )| 5 (-))+C 4 (5(+)|0(x 1 ,x 2 )| (7 (+)g(+)), (6) 

where ± in brackets indicate the helicity of partons. It should be noted that there are possible spectators. 
E.g., in the first term, the spectators must carry the total helicity A s = which can be a quark pair 
or gluon pair. In the second term the spectators must have A s = —1 because we have here A = —1/2. 
For the last two terms we have A s = ±1/2, respectively. The contributions from spectators give overall 
factors as products of ^-functions for each term. These overall factors are contained in the constants Ci 
which also depend on q in Eq.©. Because of existence of spectator-partons, the parton states in the 
above are not necessarily with the total momentum P. 

Beyond tree-level, the four matrix elements of partons in Eq.([6|) receive corrections. These corrections 
make T_ //-dependent. We define the four matrix elements as: 

T^ qg (x 1 ,x 2 ,y ,z ) = (q(p,-)\0(x 1 ,x 2 )\q(p 1 ,+)g(p 2 ,-)), 
T +qg (x 1 ,x 2 ,y ,z ) = (q(.Pi,-)g{P2,+)\0(x 1 ,x 2 )\q(p,+)}, 
T- q q(x 1 ,x 2 ,y ,z ) = {q(p 2 ,-)q(p 1 ,-)\0(x 1 ,x 2 )\g(p,-)), 

T +q q(xi,x 2 ,y ,Zo) = (g(p,+)\0{xi,x 2 )\q(pi,+)q{p 2 ,+)), (7) 

with the momenta: 

p+ = y P+ 1 p+ = (z -y )P + , p + = z P + . (8) 

The color of the state \qg) is the same as the single-quark state. Details can be found in |16} 118] . The 
four matrix elements are not independent. They are pair- wise related: 

T- qg (x 1 ,x 2 ,y ,z ) = T +qg (x 2 ,x 1 ,y ,z ), T^ qq -(x u x 2 , yo, z ) = T +qq -(x 2 , x 1 , y , z ). (9) 

As we will see, the //-dependence of the four matrix elements determine the kernel T— in different regions 
of £i i2 . Because of this, the determined kernel will not depend on the states, i.e., the coefficients Cj and 
Ci. To determine the kernel in the full region of ^i j2 , one also needs to calculate T_ of the state |n)— 
The charge-conjugated state of |n). However, T_ of \n) can be obtained from T_ of |n) through charge 
conjugation. 
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Figure 1: A set of diagrams of one- loop corrections to the denned twist-3 matrix elements. This set only 
contains the self-energy corrections and corrections from the gluon emission from a gauge link represented 
by a double line. 

We take T- qg as an example to show how the corresponding contribution to J 7 , is determined. At 
tree-level, we have: 

T^ g (x 1 ,x 2 ,y ,z ) = 2irg s y / 2z y (Nc - l){z - yo)5{z - x 1 )5(x 2 - y ). (10) 

It should be noted zq > yo because of p\ > 0. If we have the one-loop result T^j g , we can substitute the 
results of into Eq.(j3|) through Eq.©. Then we can find for £2 > 0: 

1 d 

The ^-function #(6 — £2) appears because of zq — yo > 0. Therefore, T- in the region of £1 > £2 and 
£2 > is determined by Ti qg . Similarly, one can find that in the region of £2 > 6 and 6 > is 

determined by . Combining the two contributions one has the kernel in the region of positive 6,2- 

The one-loop correction to T- qg are from the diagrams given in Fig.l and Fig. 2. The calculation 
of these diagrams is rather standard. Therefore we give the result directly. We introduce the following 
function which is just the kernel in the region of £12 > 0: 



-^1(^1,^2,6,6) = 0(x 2 ) 



£(6 ~ xi)h{£ 2 ,xi,x 2 ) +5(6 -Z2-X1 +x 2 )f (^ 2 ,xi,x 2 ) 



+<5(6 - x 2 )f 2 (Cl,xi,x 2 ] 



+ (xi o x 2 ,6 6), (12) 
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Figure 2: Another set of diagrams for one- loop corrections of the defined twist-3 matrix elements. 



with 



f 2 (£,X!,X 2 ) 



1 x 2 6(Z-x 2 ) 
2N C i{x 2 - £)+ 



9(xi - x 2 )9{x 2 - - 6{x 2 - xx)9{£ - x 2 ] 



x\ - x 2 



1 



2 £ V(£-x 2 ) 

N 2 - 1 

In x 2 — N c In I x\ — x 2 1 



(x 2 - £)+ (xi - £) 2 x x - £ 
3(N 2 - 1) 



4N r 



4(£ - x 2 ) 



0(si - 0^(^2 - si) X ) t X \ (Sx 2 - 2xi - 20 



Xl 



-20(£ - xi)0(x 2 - xi)-^(x 2 - xi - 

KX2 



0(t-x 1 )0(x 1 -x 2 ) 
1 



X2^ ~ X 2 

2x 2 xi - 3£x 2 - 2x\ + xi£ + 2£ 2 



AT2 

c 



2iV„ 



-0(x 2 - xi)6»(xi) 



-0(xi) 



Xj(x 2 - Xj) 

x§(£-x 2 ) 

^ - x 2 + Xl 



" z 2 ) 



0(-Xl)0(£-X 2 + Xl) 



X 2 {i - X 2 

The +-distributions here are defined as: 



xi#(x 2 -£-xi) x 2 -xi 

H ; 0(f -x 2 + xi. 



!i (x-z). 

(-2-X). 



?^2 
/(*) " /(*) 



+ /(x) lnx, 



x - z 

1 dz f{z) ~ f{x) +/(x)m(l-x). 



(13) 



(14) 

/a; # X 

It should be noted that our +-distribution is not the standard +-distribution defined later. The function 
/o and fi are determined by the contributions from Fig.l. In general for each diagram with a gluon 
emitted from the gauge link there is a light-cone singularity. But, the sum is free from the singularity. 
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Now we turn to the contributions from T_ qq . At tree-level we have 



r|° ) g(xi,x 2 ,?/o,2;o) = -2ng s (Nc - l)^2y (z - y )z S(x 1 - y )5(x 2 + z - y ). (15) 
With the one-loop result of T_ qq one can determine the kernel 



3 s a s ^_(xi,x 2 ,6,6)#(£i)#(-6) 



2(JV c 2-l)(£ 1 -£ 2 V=2££01n/x -« 



T^hx 1 ,x 2 ,£i,£i - 6), (16) 



The kernel determined from T_ g g is in the region of 6 > and 6 < 0. From T +qq the kernel in the 
region £1 < and £2 > can be obtained. The one-loop correction to T_ qq is still given by the diagrams 
in Fig.l and Fig. 2, where the quark line in the left side should be taken as an outgoing anti-quark. The 
calculation of these diagrams are again straightforward. The result can be summarized by the function 

Fi(xx, X2,£i, 6) = <K£l - xi)9(x 1 )(hi(—^2,xi,x 2 ) + h 2 (—^2,-X2,-xi)j 
+5(^2 - x 2 )9(-x 2 ) (hx{ix, -x 2 , -xx) + h 2 (£x, xx,x 2 )j 

-&-xx + x 2 ) yO(-x 2 )h (-£ 2 , xx, x 2 ) + 9(xx)h (£x, ~x 2 , -xx)j , (17) 



with the functions 

h (€,xx,x 2 ) = 
hx(€,xx,x 2 ) = 



2N t 



1 e(x 2 + o- X2 



0(-x 2 )6(x 2 + O 



£ + 2xx - x 2 
(£ + xi) 2 



2x9 



^ + Xx)J (£ + X 2 



xx - x 2 



+6(xx - x 2 )6(x 2 ) 



1 



+-9(-x 2 )5(x 2 + 



£ + 2xi - x 2 xx - x 2 

ic + xi) 2 i + x 2 

3 iV 2 - 1 iV 2 - 
4 



1 



h 2 (£,xx,x 2 ) 



Nr. 



1 



+ 



4 (£-x 2 ) 2 

(-xx)0(xx - 

1 £ + xx - x 2 
2N C (£-x 2 ) 2 



In |x2| — iV c In |xi — X2I 
2£ 2 - xxi + 2x\ + 3£x 2 - 2xix 2 



+0(-xi)0(zi - x 2 )- — — (2£ + 2xi - 3x 2 ) 



X-2 



+ - 



Ni 



2N r . 



-9(-x\)9(xx - x 2 ) 



9{-xx)9{i-x 2 + xx) 
xx{xx - x 2 ) 



at \ Xl . at , i~x 2 + xx 
9[xx - x 2 ) V 9\X2 - xx) 

X2 £ 



x 2 (C-x 2 ) 



(18) 



The function J- 2 is just the kernel in the region of £1 > and £2 < 0. The kernel in the region of £1 < 
and £2 > can be obtained from T +qq through the relation in Eq.Q. 
Finally, we have the kernel in the full region of £1,2 as: 



T-(xx,X2,ti,€ 2 ) 



9(tix)9(-&)Mxux 2 ,t 1 ,t 2 ) + 0(6)0(-£i) Jbfo, xx, 6, 6) 
+0(6)0(6)^1, *2, £1,6) + 9(-tx)9(-&)J r i (-X2, -x x ,-i2, -&). (19) 
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The kernel in the region with £12 < can be obtained from that in the region with £^2 > through 
charge-conjugation. Eq. (j!9[) is our main result. Since Tp and Ta are not mixed under renormalization, 
the correct kernel should satifiy: 

J 7 _(xi,x 2 ,6 1 6) = .F-0c2,£i,6,£O- (20) 

Our result has the property. Therefore, taking the part of T- symmetric between x\ and X2, one can 
obtain the evolution of Tp as a convolution with Tp. The anti-symmetric part of T"_ gives the kernel of 
the evolution of Ta. 

The evolution kernel for arbitrary x\ 2 is rather lengthy. But the convolution in Eq.(|4]) is a one- 
dimensional integral at one-loop. In some special cases, the kernel is simplified. E.g., for x\ = X2 = x > 0, 
i.e., the case where the gluon carries zero momentum fraction entering a hard scattering, we have: 

d ~ , , a.f + . f 1 dz 1 



d In fi 



Tp( W ) = 7{ - ^(^,.) + / ^ (jV c f F (x,0 - ^(«) 



d - 

T A {x,x,fi) = (21) 



dln/j 

with z = x/£. Here the +-distribution is the standard one defined as: 

1 0(z-x) , s f 1 tiz) - til) . . . . 
dzji M{z)= dz—^ — +t(l)ln(l -x). (22) 

(1 _ z )+ Jx 1 — Z 

The result in Eq. (|21|) agrees with that in [8j . There are also cases in which a quark carries zero momentum 
fraction entering a hard scattering. In these cases, one has either x\ = or x<i = 0. The evolutions in 
these cases are: 



d nn , n v a s ( f 1 dz 
-T-ip,x,n) 



din fi it [ Jx z 



1 / N ~ 1 ~ 

+ (1^ K-(0,{) + f T_(x -{,«)- — TAi - 

Changing every T_(xi,X2) to T_(x2,xi) in the above equation, we obtain the evolution of T_(x,0). 

The two studied twist-3 distributions do not mix in the evolution with other twist-3 distributions 
mentioned in the beginning, e.g., with e(x) and hi(x). But the evolution of the distribution e(x) or 
hi{x) does mix with Tp and Ta[I2J [131 El- O ne ma y derive the evolution of e(x) and hi(x) through 
their relations to Tp and T\, respectively. But the derivation is very tedious. It has been shown in |12j 
that in the large- N c limit the evolution of e(x) or Hl(x) is simplified and obeys DGLAP-type equations. 
The result for this in momentum fraction space can be found in |14j . These evolutions in the large- iV c 
limit can also be derived with the approach used here. We take e(x) as an example to give our result. 
e[x) is defined in the light-cone gauge as: 

eix) = P + I — e- iyxP+ (P\4>(yn)4>(0)\P). (24) 

J 47T 
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To derive the evolution of the two-parton twist-3 distribution in the large- iV c limit one can use the same 
approach for deriving the evolution of three-parton twist-3 distributions. The derivation is straightfor- 
ward. We have for x > 0: 



de(x) 
din fx 



a s N c 
2vr 



(25) 



with z = x/£. This result agrees with that in |14j . Although the evolution of e(x) and Hl(x) gets 
simplified in the large- N c limit, the evolutions of Tp,A are n °t much simplified in the limit. This can be 
observed from our results. This supports the observation made in [19] that simplifications of evolutions 
of higher-twist distributions can be accidental. 

To summarize: We have calculated the twist-3 distributions of multi-parton states. The distributions 
are defined with chirality-odd quark-gluon operators. The evolution kernel of the distributions has been 
obtained by the calculation at one-loop. The evolution at one-loop is an one-dimensional convolution. 
In some special cases the evolution takes a short form. We have also derived the evolution of e(x) in the 
large- N c limit and found an agreement with existing results. 
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